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1 Introduction 



Self-similar models, with a group acting transitively on 3-dimensional hypersurfaces, has been 
discussed in the literature for several decades. Within this family of models one finds a number 
of different interesting physical phenomena; chock waves and violation of cosmic censorship being 
perhaps the most prominent ones. Self-similar models are also interesting because they constitute 
asymptotic states for more general models and that they thus act as building blocks when it comes 
to understanding wider classes of models. They also generalize the spatially homogeneous (SH) 
models which admit two spacelike commuting Killing vectors, and in this context they can be seen 
as the first step towards the construction of inhomogeneous cosmological models. They might also 
shed light on the generality of properties found in the SH subclass of models. For example, in this 
article we discuss if whimper singularities, found in some of the SH models, exist in the self-similar 
case. 

One finds two approaches toward the self-similar perfect fluid models. One is the "fluid adapted" 
approach in which one choses a timelike coordinate along the fluid lines (see e.g., Q-Q). The second 
is the "homothetic" approach in which one chooses a coordinate along the orbits of the homothetic 
Killing vector (see e.g., |^]~|^]). In general, the symmetry surface changes causality. 

One often considers models that are diagonalizable. In this case the homothetic approach, taken 
together with the diagonal gauge requirement, has the drawback that one has to cover the spacetime 
with several coordinate patches. On the other hand, the homothetic approach reveals that there 
is a considerable structural similarity between the field equations for the self-similar models and 
the hypersurface-homogeneous models (e.g., SH models). There is a considerable literature about 
on how one deals with hypersurface-homogeneous models (much larger than the one on self-similar 
models). Thus the homothetic approach makes it possible to transfer ideas from the hypersurface- 
homogeneous arena to the self-similar one. Because of this advantage we will use the homothetic 
approach. However, note that results obtained in this picture can be transferred to the fluid adapted 
formulation by a coordinate transformation and vice versa. 

In this article we will focus on the "spatial part" of self-similar models exhibiting a locally 
rotationally symmetric (LRS) isometry group (in addition to the self-similar symmetry). The line 
elements for the spatially-self-similar (SSS) LRS models have been given by Wu [^. Collectively 
they can be written as 



-dt^ + Di{tfdx^ + D2{tfe-^^'= (dy"^ + k'^ sin {y/ky)dz^y^ (1) 



where /, a, k are parameters describing the symmetry groups of the various models. Canonical 
values for these parameters are given in Table 1. Note that o/c = and that, for the sake of brevity, 
we have denoted the LRS type fl (jV) models with fl (/V) even though they also contain a type 
/VIIq (/VII/i) group (see @). The SSS models with spherical symmetry are denoted by *KS where 
KS stands for Kantowski-Sachs in analogy with the corresponding SH case. 

We will consider perfect fluid models. The energy momentum tensor. Tab, is thus given by 

Tab = jj-UaUb + p{UaUb + Qab) , (2) 

where jl is the energy density; p is the pressure; and the 4-velocity of the fluid. We will assume 

p = {j-l)ll (3) 

as an equation of state where the parameter 7 takes values in the interval 1 < 7 < 2, which includes 
dust (7 = 1) and radiation (7 = 4/3). Thus we have excluded the value 7 = 2, which corresponds to 
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Table 1: Canonical choices of the symmetry parameters a, f and k. 



a stiff fluid. The reason for this is that the corresponding models behave quite differently compared 
to those in the interval 1 < 7 < 2, and thus need special treatment. Note that all dust models are 
known [10|. Note also that all self-similar LRS models are of Petrov type D (or 0) and that the 
magnetic part of the Weyl tensor is zero. 

The outline of the article is the following: In section 2 we rewrite the field equations in two steps. 
We first express the field equations in a set of variables associated with the normal congruence of 
the symmetry surface. We then introduce a dimensionless set of variables in order to obtain a 
maximal reduction of the coupled system of ordinary differential equations. 

In section 3 the reduced phase spaces and invariant submanifolds of the various models are 
discussed. The dimensionality of the fully reduced phase spaces together with the relation between 
the models, in terms of Lie contractions, are given in a diagram. Similarities and differences between 
the various models are discussed. The relation between the field equations for the SSS models and 
the timelike-self-similar (TSS) models and the possibility of extending the SSS models to the TSS 
sector is commented on. 

In section 4 an equilibrium (critical, singular) point analysis of the SSS type il and LRS 
models is carried out. Asymptotic expressions for the line element and the kinematic fiuid quantities 
are given. The remaining SSS and TSS LRS models need special treatment, and will be discussed 
elsewhere. 

In section 5 complete phase portraits are given for a number of type il and jY models. Some 
global issues are also discussed. In section 6 we discuss the nature of the singularities occurring in 
the SSS type il and jV models. Appendix A describes the properties of the fiuid congruence and 
the condition for the spacetimes to belong to Petrov type 0. Appendix B gives the relation between 
the presently used coordinates and those of the fiuid approach. 



2 Derivation of the autonomous DE 

Misner has introduced a useful metric parametrization in the context of SH cosmology [11 1. In the 
present case an analogous parametrization can be introduced 

D, = e^'-^^* , D2 = e^'^"^ . (4) 

This parametrization is closely related to the kinematic properties of the normal congruence of the 
ds^ geometry (which is conformally related to the SSS geometry). The expansion, 9, and the shear. 
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which can be described by a quantity cr+, of this congruence, are related to /J^'"*" by 

= 3/3°, a+ = 3/3+ , (5) 

where the dot stands for d/dt. The tetrad components of the fluid velocity are conveniently 
parametrized by 0, 0)/Vl — v"^ where v is just the 3- velocity with respect to the symmetry 
surfaces. 

Einstein's equations. Gab = Tab, and T^^.^ = lead to: 



Evolution equations 
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1 



3 3 
(7+ = -ea+ + 2 



ai + 2ia + f)fBt-- 



1 (37 - 2) + (2 - 7) 



2 1 + (7 - 1) 



vBi + 2{a + f)fBl-kBi , 



Bi = -{-e + 2a+)Bi , 



B2 = --{e + a+)B2 , 



1 -t;^ 
37(1 -(7-1)^^2) 



' ((37 - 4) + 2(7+) ^; + 3 ((7 - 1) (2/ - (2a + 3/) 7) v'' + {2 - 7) /) i?i 



where 



Bi = e 



£>2 



(6) 
(7) 



The quantity ^„ stands for the energy density of the fluid, measured by an observer associated with 
the normal congruence of the symmetry surfaces, multiplied with the factor e~^-^^. It is related to 
fl, the energy density of the fluid, by 



1 + (7 - ly _ 



-fie 



2fx 



The above evolution equations are not independent. Instead they are related by 

Constraint equation 



(8) 



-fviin + 2 l + (7-l)i;^ 



a + =-f]a+-l-fe 



Bi = . 



Defining equation for /in 

fin = l{e^-(Tl-S (3a + /) (a + /) Bj + 3kBl) . 



(9) 



(10) 



Of interest is also: 



Auxiliary equation 
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Z^" = i + (^_i)^2 (^-^ + 2(« + f)Biv + - {2a+ - e) v'j . (11) 

Note the close relationship between the above presentation of the field equations and the one 
given by Hewitt and Wainwright for the SH type V case |12]. This close connection will allow us 
to transfer some of the ideas used by Hewitt and Wainwright to the present class of models. In 
particular, we will follow their treatment of the constraint and the introduction of ^-normalized 
dimensionless variables: 

The density Hn is replaced by the density parameter which is defined by 

f^n = ^ . (13) 



An introduction of a dimensionless time variable r, 



dt _3 



(14) 



leads to a decoupling of the ^-equation 



flf) 

^ = -{l + q)9, (15) 

where the deceleration parameter q, associated with the normal congruence of the SH ds'^ geometry, 
is given by 

g = 2S+ - 6(a + f)fA + ^ ^ _ -^^^2 ^« • (16) 

Alternatively one can write q as 



,2,..^ l(3(2-7) + (57-6) 



2\ 



g = 2 - 6 (a + /) (3a + 2/) A' + 6kK - ' '\ ' ' ' ' , (17) 

2 1 + (7 — Ij 

by using the definition of r2„ . The remaining evolution equations can now be written in dimension- 
less form: 

Evolution equations 

= [-2 + g + 2vA (3a + 2/)] S+ - 2fvA + 6fA^ (a + f) - 3kK , 
A' = {2^+ + q)A, 
K' = 2{-Y.+ + q)K , 

7(l-(7-l)i;2) 



V = — 7^ [7 ((37 - 4) + 2S+] ^; + 3 ((7 - 1) (2/ - (2a + 3/) 7) v' + (2 - 7) f) A\ , 

(18) 

where a prime denotes d/dr. 

Constraint equation 
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-tv^n + 2 [l + (7 - 1) [(3a + 2/) S+ - /] A = . (19) 
Defining equation for ri„ 
Jl„ = l-S^-3(3a + /)(a + /)^2^3A;K . (20) 
Auxiliary equation 



0' 



-(37 - 2) + 2g + 67(0 + f)vA + (27S+ + 2g(7 - 1) - (2 - 7)) V 



(21) 



- 1 + (7 - l)t,2 

The above equation system exhibits a discrete symmetry. The field equations are invariant 
under the transformations 

(S+, ^, K, v) ^ (S+, K, -v) , . (22) 
Tlie line element can be obtained when A and 9 has been found through the relations 

Dl = {OA)-^ , Dl = (e^K)-' , i = 3/ y . (23) 

The relationship between the kinematic fluid quantities and the 9, K, A,Ti^,v variables is given in 
Appendix A. 

3 Reduced phase spaces and invariant submanifold structure 

The K-equation decouples for the type il and jV models leaving a reduced system of equations for 
the u variables, related by the constraint in eq.(^). The relationship between the various 

models is given in terms of Lie contractions in Diagram 1 together with the dimension of the 
reduced phase space. 

We have chosen to present the equation in expansion normalized variables since this leads to 
relatively simple equations. However, a very desirable property is compactness of the reduced phase 
space. In this context eq.([lO|) plays an essential role. One can use the fact that fin is nonnegative 
to produce inequalities which can be used in order to find "dominant" quantities. The expansion, 
9, is such a quantity for some of the models (those for which k < 0, (3a + /)(a + /) > 0; as seen 
by eq.(|To|)). Note that 9 cannot change sign for the matter filled "0-dominated" models because of 
eq. ([To|) . Furthermore, because of the discrete symmetry in eq.([2^) one can therefore, without loss 
of generality, assume A > (this just corresponds to looking at models with positive expansion). 
Unfortunately one cannot make this assumption for the other models where 9, and thus also A, 
may change sign. 

The boundary consists of a number of invariant subsets which plays an important role when 
determining the qualitative properties of the orbits in the reduced phase space. The boundary of 
the reduced phase space will thus be included. We have the following invariant submanifolds on 
the boundary: (i) A = 0, (ii) v = 1, (iii) v = —1, (iv) the vacuum submanifold 1 — — 3(3a + 
/)(a + f )A'^ = 0. For the k models we also have (v) K = as an invariant submanifold of the 
reduced phase space. 

The V = ±1 submanifolds has a physical interpretation in terms of models with directed fluxes 
of neutrinos (see [^). In the vacuum submanifold case one can interpret v as the velocity of a test 
fluid (see [^). The constraint leads to vCln = when A = 0. For v = we obtain the reduced 
equations for the orthogonal SH LRS type I, III and KS models. When fi^ = one obtains the 
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Diagram 1: Lie contractions for self similar and spatially homogeneous LRS models. The abbre- 
viations Cans., HH, HSS stand for causal character, hypersurface homogeneous and hypersurface 
self-similar respectively. 



reduced equations for the vacuum models for the same models but with a test fluid. For the k ^ 
models, the K = manifold yields the same equations as the reduced type il equations. Note that 
f = 0, in contrast to the SH models, is not an invariant submanifold (except for 7 = 2). Thus 
one naturally obtains, for the present class of LRS models, the well known result that SSS models 
are tilted if 7 / 2 (i.e., there are no SSS models with a fluid flow orthogonal with respect to the 
symmetry surfaces) p. 

3.1 The reduced phase space of type jV and il models 

For these models k = 0. This results in a decoupling of the equation for K. Since K decouples the 
boundary value ii' = is of no interest (it is enough to analyze the reduced system of equations). 
Thus K satisfy the inequality K > 0. This leaves a coupled system of equations for the remaining 
variables (S+, A, v), related by the constraint equation. As mentioned above, the (3a+/)(a+/) = 
and the (3a + /)(a + /) < models, do not have a compact reduced phase space. Moreover, the 
3a + / = and a + / = models behaves drastically different compared to the other models. 
This is seen in the equilibrium point analysis below. Shikin has been able to find the general 
solution for the 3a + / and a + / = cases [Q, |[. Shikin has also qualitatively investigated the 
(3a + /)(a + /) < models [|. Below we will focus on the "6'-dominant" (3a + /)(a + /) > 
case which contains the SH type V models and the SSS type il models. Recall that one can 
assume that ^ > for the ^-dominant models. However, we will also list the equilibrium points 
for the remaining models since this gives an indication of how different the three types of models 
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((3a + /)(a + /) < 0, (3a + /)(« + /) = 0, (3a + /)(a + /) > 0) are. 

There are four invariant submanifolds describing the boundary of the reduced equation system 
of the ^-dominant models: The v = ±1 submanifolds; the vacuum submanifold 1 — — 3(3a + 
/)(a + f)A^ = 0, and the ^ = submanifold which can be divided into three parts depending on 
u < 0, f = 0, w > 0. 



3.2 The reduced phase space of type ^III and *KS models 

It follows from eq.([lO|) that one obtains a compact reduced phase space in the (SSS) type models. 
However, the 9 normalized variables are not compact for the *KS models. In a forthcoming article 
we will show how one can use eq.(|To|) to produce compact reduced phase space variables in the *KS 
case. The type JIII models are characterized by the choices k = — l,a = 0,/ = —1. The reduced 
system consists of equations for the T,j^^A,K,v variables, related by the constraint in eq. ([l9|). 
Here the K = {) value is included. This leads to the inequality K >Q. The boundary is described 
by five pieces: The v = ±1 submanifolds; the vacuum submanifold 1 — — 3^4^ — ?>K = 0; the 
^ = submanifold, and the K = Q submanifold. 



3.3 The relation between the SSS and TSS field equations 

It is important to note that the above equation system (where /i„ (r2„) is assumed to be eliminated 
by use of eq.([To|) (eq. (pO|) )) also describes the field equations in the TSS region if one makes the 
change K —K in the k ^ case. However, in the TSS case > 1 instead of, as presently, 
< 1. 

It's natural to extend solutions to the TSS sector when the symmetry surface changes causality 
and when one does not run into a (non-coordinate) singularity (compare with the LRS type V 
discussion by Collins and Ellis ||l^). A number of phenomena occurs in the TSS region which do 
not exist in the SSS case. Hence, we will investigate the TSS region separately in a forthcoming 
article. This will allow us to obtain a more global picture of the self-similar LRS models than in 
the present article which will, from now on, deal exclusively with the SSS case. 



4 Equilibrium points and asymptotic analysis 

The reduced phase space is determined by a coupled system dX/dr = F(X) subject to a constraint 
G(X) = given by eq.([l9|) (X constitute the reduced phase space variables). Of central importance 
to the investigation of the dynamical system are the equilibrium points which are determined by 
the equations F(X) = 0, G(X) = 0. Below the equilibrium points are given together with which 
shows if an equilibrium point is located on a vacuum submanifold or not. The gradient of G'(X), 
which is used to locally solve the constraint to linear order, and the eigenvalues and eigenvectors 
of the remaining locally unconstrained system are then listed. 

A local analysis of the equilibrium points is then used to express the fluid energy density (/i), 
expansion (6), shear scalar (a), acceleration (a), divergence of the acceleration (a^-^); the line 
element (Di^2); and the Weyl scalar (C), for the various asymptotes in the interior matter filled 
phase space. The tables describing asymptotic properties differ from the ones given by Collins and 
Ellis for the SH type V models [13|. Their tables relate the initial and final state of a given orbit (or 



orbits). Such tables will be given in the next section when we discuss global behaviour. Our tables 
also differ from those in reference [O] in that we express the behaviour in terms of the independent 
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variable r instead of the proper time associated with the comoving fluid frame. The reason for this 
is that the SSS and SH models need separate treatment if one wants to use proper fluid time. The 
results in this section and those by |13, ^] for SH type V models are easily related by using, e.g., 
6, as independent variable. We only give the r-dependence; the self-similar dependence on e-^^ is 
given in Appendix A. The acceleration scalar, a, and the divergence of the acceleration, a|^, are of 
course zero for the dust, 7 = 1, models. 



4.1 Type jV and il models 

The reduced phase space consists of the constraint surface in the A, i;) space, subject to the 
inequalities Qn > 0,u^ < 1. For the ^-dominated models we also have A>0. 



4.1.1 Equilibrium points with zero tilt, v = 



The equilibrium point F: 

S+ = 0, ^ = 0; 0„ = 1; 

VG = (0, -2/, 7) , {v eliminated) ; (24) 

-1(2-7), (1,0); i(37-2), (0,1), 

for 1 < 7 < 2, where the constraint has been used to solve for the v variable and where the 
eigenvalues of the resulting unconstrained system are grouped together with the corresponding 
eigenvectors. The equilibrium point is a saddle point. In Table 2 we give asymptotic expressions 
for the single initial asymptote entering the interior of the reduced physical phase space of the type 
jV and il models. 
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Table 2: Asymptotic expressions for the single initial asymptote coming from the equilibrium point 
F. Here, and below As. St. is an abbreviation for asymptotic state. Initial states are denoted by i 
and final states by f. 

The equilibrium points Kj_: 

S+ = ±1 , A = 0; nn = ; 

VG = (0,2(3a + /) ,0) , for S+ = +1 , 

(25) 

VG = (0, -6 (a + /) , 0) , for S+ = -1 ; (^eliminated) ; 

3(2-7), (1,0); 2S+ + (37-4), (0,1), 

for 1 < 7 < 2. Different equilibrium points are often closely related through sign changes or zero 
values of v and In such cases we denote the equilibrium points with a common kernel and 
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indicate the value of v in the upper position and of in the lower (if there only exists points with 
a given value of v or S+ we do not indicate the corresponding value). Thus in the present case of 
K^, refers to f = and it refers to the sign of (Note that this leads to a somewhat different 
notation than that used by Hewitt and Wainwright for the SH type V case The point is a 

saddle point on the boundary and does not give rise to any asymptotes in the interior phase space. 
The point on the other hand is a stable source giving rise to a 1-parameter set of asymptotes 
described in Table 3. 
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Table 3: Asymptotic expressions for the fluid congruence around the i^9-point. 



The equilibrium point M^: 

This point only exists for / = and therefore we can choose A > 0. 







1 

3a 



; 



VG = (2, 0, 0) ; (S+ eliminated) ; 

-(37-2), (1,0); 37-4, (0,1) 



(26) 



There exists a bifurcation for 7 = 4/3, which needs special consideration; we refer to [|l^]. For 
1 < 7 < 4/3 the point is stable leading to a 1-parameter set of final asymptotes, among them 
the Friedmann asymptote. For 4/3 < 7 < 2 the point is a saddle and has no asymptotes into the 
physical phase space, except for the special Friedmann asymptote. The asymptotic behaviour is 
described in Table 4. 
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Table 4: Asymptotic expressions for the fluid congruence around the M^'-point. 
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4.1.2 Points with intermediate tilt, < < 1 



The equilibrium points M^±: 1< 7 < 2; 7 / 2//(2a + 3/). 



V / , _ e „, _ nc2 ^ / 7^ci + 4(2-7)(7-l)/ci ^ o -n- 

+ "3a + 2/' 3a + 2/' "^"201(7-!) V 4c2(7 - 1)' ' 

ci := (37 - 2)/ + 2a7 , C2 := 2(7 - 1)/ + (37 - 4)a ; 

/ 2ei;^(3a + 2/) (7 - 1) - 2/7^ + 2e(3a + 2/) Qe^v{3a + f ){a + f ) \ 
\^ 3a + 2/ ' 3a + 2/ ' ) ' 

(^eliminated) ; - " ^ "n " "^^^ ^ ' 

(3a + 2/)7v (7 - Ij^cf 

ei := -3e((2(7 - 1)(2 - ^ff - 4f^'){^ - l)c? + 2(4(7 - 1)' - (2 - 7))(7 - l)'c?/ 

+ 4(7 - rfa-f4 - 4(7 - 1)(7 - 2)cici/72 + ch^) , 
62 := -3((3(7 - 2)/ + ci7)(7 - 1)(7 - 2)ci/ - (7 - ifcf - (7 - 2)4fj^)c27 • 

(27) 

The dr in ikf* refers to the different signs of the root in the expression for v. The condition 
7 7^ 2//(2a + 3/) corresponds to ci 7^ 0. The eigenvector of the second eigenvalue above points in 
the direction of the vacuum submanifold, while the eigenvector of the first eigenvalue points into 
the physical phase space. The expression for the eigenvector associated with the first eigenvalue 
is quite cumbersome and will not be given. Note that there are two sets of points for the non-0- 
dominated models since A can be both positive and negative in that case. For the ^-dominated 
models there is only a single pair of points. The requirement of A > taken together with the sign 
of 3a + / determines the sign of the parameter e which can take the values ±1. The first eigenvalue 
is negative for the ^-dominated models. Hence it follows that, in the ^-dominated cases, possible 
initial asymptotes lies in the vacuum submanifold. Thus the only interior asymptote is a final one. 
Its generic behaviour is given in Table 5. 
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Table 5: Asymptotic expressions for the fluid congruence around the M^^^ -points. The constants 
ri,2 are given by Fi = g^'^^/ and F2 = ^^a+2f ■ Here and below will stand for a constant which 
will not be given because it is quite complicated (the actual value of is different for different 
tables) . 

The equilibrium points M^±: 7=1. 
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f If 

A = ± -, v = ±-; n„ = 0; 



3a + 2/ ' 3a + 2/ ' a 

^ / 2(3a-/) (a + /) 6/ (3a + /) (a + /) \ (28) 
V a(3a + 2/) ' a (3a + 2/) '7 ' 

(A eliminated) , 

The local analysis yields a single multiple valued eigenvalue (m below) and a linear solution of the 
following form: 



^ 3a + 2/ " ' a 

(29) 

3(a + /)(a-/) _ 3/(3a + 2/)(a-/)(a + /)^ 
(3a + 2/)a ' (3a - /)a3 

where /ci and k2 are constants. For ^-dominated models with (3a + /) < and models with 
(a + /) > and a — / < none of these points exist in the physical phase space. However, for 
models with (a + /) > and a — f > the point M'"+ exists and is a sink. For non-0-dominated 
models it can be a sink or a source depending on the values of a and /. Asymptotic behaviour is 
given in Table 6. 
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Table 6: Asymptotic expressions for the fluid congruence around the -points for 7 = 1. The 
constants Fi^a are given by Fi = 3^ and Fa = while Af = 

T/ie equilibrium points M^*: 7 = 2//(2a + 3/). 

- - ^ A = ± . - = ±^^; ^n = 0; 



^ 3a + 2/ ' 3a + 2/ ' 4a + / 

^ / , 4(6a^ + 2a/-3/2) (3a + /) 12/(3a + /)(a + /) \ 

1,^ (4a + /)2(3a + 2/) ' (4a + /)(3a + 2/) ' ' ^^"^ 

(^eliminated); 12(a + /)(a-/) . 6(4a + /)(a + /)(a - /) 



(3a + 2/)(2a + 3/) ' (Sa^ + 4a/ + /2) (2a + 3/) ' 

This point only exists for models which are non-0-dominated, and therefore wc have only given 
the eigenvalues and not the eigenvectors. We have also refrained from giving a table with the 
asymptotic properties. 
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4.1.3 Points with extreme tilt, 

The equilibrium points M^: 

/ 



3a + 2/ ' 



A = T 



3a + 2/ 



±1 ; nn = 0; 



3a + 2/ ' 
6(a + /)(57 



3a + 2/ 

12(a + /) 
(3a + 2/) ' 



(S+ eliminated) ; 



(31) 



(1,0) ; 



6) 



(3a + 2/)(2-7) 



(0,1) 



For ^-dominated models with (3a + /) < only M"*" exists, while for models with (a + /) > 
only M~ exists. The quotient (a + /)/(3a + 2/), appearing in both eigenvalues, is positive for all 
6'-dominated models. However, for non-^-dominated models different signs are possible and thus 
this factor will affect the stability of the point in this case. For ^-dominated models, the only 
bifurcation is for 7 = 6/5. For 7 < 6/5 the existing point is a saddle and there is no asymptote 
from or into the interior region of the phase space, while for 7 > 6/5, the points are always sinks. 
Note that the fluid quantities require separate treatment when 7 = 4/3, even though there is no 
bifurcation (note the asymptotic behaviour of a in Table 7). 

For 7 > 6/5, the sink associated with the existing point does not correspond to a curvature 
singularity. For ^-dominated models with (3a + /) < there is no final singular behaviour at all. 
This is also true for the case (a + /) > and / < — 2a(37 — 4)/(57 — 6); but if / < — 2a(37 — 
4)/(57 — 6), then the point corresponds to a crushing singularity (i.e., the kinematic quantities blow 
up while the curvature stays finite (it goes to zero)). 

In the above discussion about singularities, we should really only say that the curvature and/or 
kinematic quantities blow up asymptotically since we have not investigated if one really reaches 
the asymptotic stage in a finite amount of proper fluid time or in a finite affine distance along a 
geodesic. Such an investigation will not be undertaken in the initial and final asymptotic analysis 
below either. To prove that we really have singularities we would need to do an additional analysis 
involving asymptotic coordinate variable transformations for each case. However, for the closely 
related SH case and for the unphysical SH geometry associated with the SSS geometry, the initial 
and final states corresponding to possible singularities are reached in a finite time. Hence we believe 
that these states correspond to singularities and we will thus refer to them as such. 
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Table 7: Asymptotic expressions for the fluid congruence around the M^-point. The constants 

Fi o ^ are triven bv V^ - 3((°+/)(57-6)-(2-7)a) p _ 6((37-2)/+2a7) , p _ 2a+f 
i 1,2,3 are given Oy i 1 - (3a+2/)(2-7) > ^ 2 - (3a+2/)(2-7) ^ 3 - 3a+2/- 
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The equilibrium points Ti^: 

S+ = -l±3(a + /)^ , ^ = ±1; n„ = 
27(3a + 2/ - /E+; 



2(l + S+)(/-(3a + 2/)E+) 



,±27 (3a + 2/-/S+), 



3(a + /) 

2(7-2)(l + E+)((3a + 2/)S+-/) 
3(a + /) 



3(a + /) 
(yl eliminated) ; , (1,0) ; 

-2(2(2a + 3/)S+-a) ^ ^ ^^^^ _^ _ /) e^, 3 ((2a + 3/)S+ - a) 62) , 

o + / 

ei := (2/(37 - 5) - 3a(2 - 7)) - (3a(37 - 4) + 2/(37 " 5)) S+ + 37a - (37 - 2)/ , 
62 := 3a + 2/ - /E+ . 

(32) 

For ^-dominated models only one of the lines exists in the physical part of the phase space. For 
models with (3a + /) < 0, only exists, while for models with (a+f) > only 7^+ is present. The 
condition that J7„ > 0, limits the possible values of to lie in the interval —1 < < //(3a + 2/). 
We also note that for the models with (3a + /) < the existing line of equilibrium points splits 
into two parts. The part where 11+ attain the values —1 < S+ < 2a/(2a + 3/) constitutes a source 
while the part where 2a/ {2a + 3/) < S+ < //(3a + 2/) is a sink. For models with (a + f) > 
and (f — a) < 0, the entire line 7i~^ is a source. However, for models with f > a the line splits 
into two parts where the stability of the two parts is governed by the same inequalities as for the 
(3a + /) < case. The asymptotic behaviour is given in Table 8. 

All 0-dominated models with an initial point on one of the lines Tl^ have an initial crushing 
singularity. This crushing singularity is a curvature singularity if Fi in Table 8 is negative. Thus 
for models with / > there is an initial curvature singularity. 

The part of H"^ corresponding to a sink is never associated with a final curvature singularity 
when it comes to 0-dominatcd models. If (3a + /) < and 2a/(2a + 3/) < S+ < (2a + /)/(4a + 5/) 
there are no final singularities at all, but for (3a + /) < and (2a + /)/(4a + 5/) < S+ < //(3a + 2/) 
there is a crushing singularity. For ^-dominated models with (a + /) > and a < f < 2a, there is 
no final singular behaviour at all. For 2a < /, H'^ again splits into two parts whose behaviour is 
governed by the same inequalities as in the "splitted" case. For non-0-dominated models the 
situation, which will not be discussed further, is more complicated. 



As. St. 




e 


a 


a 




Di 


D2 


C 


i,f 






gr^r 


gr^r 


g2r2r 


g(2E+-l)r 


e(i+s+)^ 


gFir 



Table 8: Asymptotic expressions for the fluid congruence around the 7^='=-points. The constants 

Fi o are ffiven bv - -Ml±5±) _ (5S+-i)/-2a(i-2S+) 
i 1,2 are given oy 1 1 - ^^^j,^ , i 2 - 
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The equilibrium points K^: 

S+ = ±1 , ^ = , v"^ = 1 ; 17„ = ; 

VG = {-2-fsgnv) , 27 (3a + /) , 0) , for S+ = 1 , 

VG = (-27sgn(t>), -67 (a + /) , 0) , for S+ = -1 ; (33) 

(^eliminated); 2(1 + S+) , (1,0); --^ (2S+ + (87 - 4)) , (0,1). 

2-7 

The points constitute one of the boundaries of the physical part of the lines TC^. They describe 
the boundary asymptotes which play an important role in the dynamical system analysis. This 
motivates their special treatment. The points have no asymptotes entering the interior region 
of the phase space. 



The equilibrium point Ai 



± . 

f 1 



^ 3a + 2/ ' 3a + 2/ ' 

0, (1,0); ,0,,, 

These points are at the other end (compared to K^) of the lines TC^. The points have no asymptotes 
entering the interior region of the phase space. 

5 Global SSS behaviour 

5.1 Type il models 

The type il models are defined by a = and / = — 1 (see Table 1). The model is ^-dominated 
which implies that our variables leads to a compact reduced phase space. Only the line TC~ of 
the two lines Ti^ lies in the physical phase space. The stability of Tl~ is governed by the sign of 
E+; for S+ < it is a source, and for S+ > it is a sink. Note that the requirement > 
limits the range of S+ to —1 < < 1/2. The only bifurcation appears for 7 = 6/5, when the 
equilibrium point M'"+ passes through M"*" and enters the physical phase space, stabilizing Af+. 
As 7 approaches 2, the point M^^ moves along the line M"*" — M'"+ — M~ , while v approaches 
0. Thus there are only two cases: 1 < 7 < 6/5 and 6/5 < 7 < 2. In both cases there are no 
interior equilibrium points, and hence no limit cycles. In addition, consideration of the orbits on 
the boundary leads to the conclusion that there are no heteroclinic cycles (compare with the SH 



type V case discussed in p^ ) 



5.1.1 Models with 1 < 7 < 6/5 

The different orbits can be written symbolically as Ti^ — > Tl^ , F — > Tl^ and — > Ti^. The inter- 
mediate evolution of orbits which define one-parameter families of solutions can be approximated 
by the following heteroclinic sequences: 



15 



(1) Kl^K+^M+^ H- 

(2) Kl-^F^n- 

(3) n- ^K^^F^n- . 

Prom Table 9, which describes the initial and final states of all asymptotes in the interior phase 
space (note the difference with the tables in the previous section where the tables referred to a single 
equilibrium point/line), it follows that all models begin with a curvature singularity and approach 
H~, with < S+ < 1/2, at late times (see Diagram 2a). None of these endpoints correspond to 
curvature singularities, although for S+ > 1/5 they are crushing singularities. 
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Table 9: Initial and final points for the type i/-models with 1 < 7 < 6/5. The constants ri^2 are 
given by Ti = -2(1 + S+) and V2 = (5E+ - 1). 

5.1.2 Models with 6/5 < 7 < 2 

The orbits are described hy ^ H' , F ^ , Kl ^ H- , Kl ^ M''+ and K° ^ M+. It 
follows from Table 10 that all models start with a curvature singularity. None of the end-points 
correspond to a curvature singularity. Points on 7i~ with S+ > 1/5 and correspond to crushing 
singularities. The remaining final points are not associated with any singular behaviour at all. The 
intermediate evolution of orbits can be approximated by the following heteroclinic sequences: 

(1) Kl^ K+ ^ M+ 

(2) Kl ^ M+ 

(3) Kl^F^n- 

(4) n- ^K^^F^n- . 
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Table 10: Initial and final points for the type i/-models with 6/5 < 7 < 2. The constants 
r^, // = 1 - 5, are given by Ti = -2(1 + = (5S+ - 1), = 3(57 " 6)/ (2(2 -7)), 

r4 = 3(37 — 2)r/(2 — 7) and = 3(a + /)r/(3a + 2/). Here and below fin. is an abbreviation for 
finite constant. 



5.2 Some comments on LRS spatially homogeneous type V models 



The SH LRS type V models have been extensively discussed in the literature ||T^, 13, 14, 15]. 
However, we will here make some remarks about the existence of monotonic functions. The SH 
type V models differ from the SSS type /V and il models in that there exists an exact solution 
corresponding to an invariant submanifold described by f = 0, T,-^- = (this is just the open 
Friedmann- Robertson- Walker (FRW) model). This separatrix divides the reduced 2-dimensional 
phase space of the SH type V models into two regions. For the type /V and il models one also has 
a separatrix leaving the equilibrium point F. However, at what equilibrium point this separatrix 
ends up at depends on 7 and /, and so far one has not been able to find the corresponding exact 
solution. The simple division of the phase space created by the FRW separatrix in the SH type V 
models allows one to find monotonic functions valid in the v > and f < regions respectively. 

The conditions A, fi^ > (for the interior matter part of the type V phase space) and the 
constraint, lead to < 0. Let us first consider f < 0. Then E_|_ (A) is a monotonically decreasing 
(increasing) function for 1 < 7 < 2. The variable v is monotonically decreasing when 4/3 < 7 < 2. 
If u > then V is monotonically decreasing if 1 < 7 < 4/3. Thus one can obtain considerable 
qualitative dynamical information because of the simple separatrix structure. Unfortunately this 
structure is not available for the SSS models. 
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Diagram 2: The reduced phase space of the type ±1 models with (a) 7 < 6/5 and (b) 6/5 < 7 < 2. 
5.3 Models with radiation, 7 = 4/3 

The general structure of the type jV models is complicated by the number of parameters. To avoid 
being lost in details we here therefore restrict ourselves to the physically interesting radiation case, 
7 = 4/3. For 7 = 4/3 we obtain a number of bifurcations for different values of the parameters a 
and /. 

5.3.1 Type /V models with 3a + / < 

The phase space has the same structure as for the type il models with 7 = 4/3. It follows from 
Diagram 3a that the initial and final end points of the orbits are described by H^, F 'H~, 

K\ ^ H', ^ My^ and K\ M+. From Table 11 it is seen that all models start with a 
curvature singularity and that the point M+ corresponds to a crushing singularity. If the final 
points on 7i~ are crushing singularities or not depends on which of the inequalities given in section 
4 the final value of satisfies. It is clear though, that the line has points corresponding to crushing 
singularities or no singularities at all. We have the following heteroclinic sequences: 

(1) K^^ K^^ M+ 

(2) Kl ^ M+ 

(3) Kl^F ^n- 
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(4) n- ^K^^F^H- 
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Table 11: Initial and final points for the type /F-models with 7 = 4/3 and (3a + /) < 0. The 
constants ri,2,3,4 are given by Ti = -2/(l + E+)/(a + /), = ((5S+ - 1)/- 2a(l - 2S+))/(a + /), 
^3 = and r4 = -6(4a + 3/)/(3a + 2/). 

5.3.2 Type /V models with a + / > and / < 

The orbits are described by M~ , F M~ , and M~ (see Diagram 36). It follows 

from Table 12 that the orbits starting from Ti^ do not have an initial curvature singularity. The 
other two types of orbits have an initial curvature singularity. The heteroclinic sequences are given 
by: 

(1) 7i+ ^ M- 

(2) ^ K^_^ F ^ M- 

(3) F ^ M- 

(4) K^_^ M- . 

As seen in Diagram 36, all orbits approach M~ for late times, and as can be seen from Table 12, 
this point corresponds to no singular behaviour at all. 

5.3.3 Type V models 

For / = 0, u = is an invariant submanifold, characterized by the separatrix F M°. The other 
orbits are characterized by TV' and K^, M~ (see Diagram 3c). Orbits starting from the 
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Table 12: Initial and final points for the type /^-models with 7 = 4/3 and (a + /)>,/ < 0. The 
constants ri,2,3,4 are given by Ti = -2/(1 + S+)/(a + /), = ((5S+ - 1)/ - 2a(l - 2S+))/(a + /), 
^3 = and r4 = -6(4a + 3/)/(3a + 2/). 



points i^^j*, and F have an initial curvature singularity while orbits starting from Ji'^ do not. The 
point M*^, which only exists for these models, has a zero eigenvalue for 7 = 4/3 and needs special 
treatment. We therefore refrain from giving a table of initial and final asymptotic expressions and 



instead refer to |13]. It follows from Table 7 that M~ is associated with no singular behaviour at 
all while Table 8 shows that the line W"*" corresponds to a crushing singularity. The intermediate 
evolution of orbits can be approximated by the following heteroclinic sequences: 

(1) ^ K° ^ F ^ MO 

(2) K%^F ^ M- 

(3) K%^K^_^M- . 



5.3.4 Type /V models vi^ith a + / > and < f < a 

In this case the SH type V point splits into the two points . These points move apart as / 
increases. When / approaches a, M'"+ approaches the point M"*". The initial and final stages of the 
orbits are given by W+ M''+ , F M''+ , K% M'"+ , Kl M"- and K% M" (see Diagram 
3(i). It follows from Table 13 that all orbits have an initial curvature singularity. Solutions ending 
in M~ have a final crushing singularity. The intermediate evolution of orbits can be approximated 
by the following heteroclinic sequences: 
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Table 13: Initial and final points for the type /F-models with 7 = 4/3 and (a + /) > 0, < / < a. 
The constants ri,2,3,4 are given by Ti = -2/(l + S+)/(a+/), T2 = ((5S+ - l)/-2a(l-2E+))/(a+ 
/), T3 = — ^^^2f ^^"^ ^4 ~ — 6(4a + 3/)/(3a + 2/). The constants X12 have long expressions which 
will not be given. 

5.3.5 Type models with a + / > and f > a 

The asymptotic properties of the orbits are described by — H'^, F — > 7^+, Kl — > 7^+, 
— My_ and Kl — >■ M~ (see Diagram 3e). Prom Table 14 it is seen that all models have an 
initial curvature singularity. Prom Table 14 it also follows that orbits ending in the point M~ have 
a crushing singularity. A solution ending on the line FL^ may have a final crushing singularity or no 
singularity at all. The existence of a final crushing singularity depends on which of the previously 
given inequalities the final value of S+ satisfies (see section 4). We have the following heteroclinic 
sequences: 

(1) Kl^K^^ M- 

(2) Kl M'"- M- 

(3) Kl^F^n-^ 

(4) Kl M"- ^ 7^+ 

(5) 7i+ ^ ^ F ^ . 

5.4 Some remarks 

All models starting from points with non-extreme tilt start from the points F or Kl and are 
characterized by curvature singularities. Models belonging to the class a + / > and / < 0, 
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Table 14: Initial and final points for the type /T^-models with 7 = 4/3 and (a + /) > 0, / > a. The 
constants ri,2,3,4 are given by Fi = -2/(l + S+)/(a + /), = ((5S+ - 1)/- 2a(l - 2E+))/(a + /), 
^3 = -i0f and r4 = -6(4a + 3/)/(3a + 2/). 

which initially have extreme tilt, may be extended so that they come from the TSS region. The 
remaining models on the other hand have an initial curvature singularity. A solution which ends 
at a point corresponding to extreme tilt can be extended to the TSS region since these points do 
not correspond to a curvature singularity (they only correspond to a crushing singularity or no 
singularity at all). There are no interior equilibrium points in any of the ^-dominated type jV 
models, and hence no limit cycles. Furthermore, consideration of the orbits on the boundary leads 
to the conclusion that there are no heteroclinic cycles either. 
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Diagram 3: The reduced phase space of the fV models with (a) 3a + / < 0, (b) a + / > and 
/ < 0, (c) a > and / = 0, (d) a + / > and < / < a and (e) a + / > and a > 0. 



23 



6 Discussion 



Note that we do not obtain any final conformal singularities (i.e., curvature singularities for which 
only the Weyl tensor blows up) at the Cauchy horizon since C — > whenever — >■ 0, const. This is 
in contrast to |]l^ where it is claimed that one has a conformal singularity at the horizon for the SH 
type V models with 4/3 < 7 < 2 (see the lower part of figure 7 in We only obtain a crushing 

singularity in this case and thus one should be able to extend the space time beyond the Cauchy 
horizon. This result is supported by our asymptotic analysis (done with help of the computer 
algebra systems REDUCE and the SHEEP package CLASSI) as well as numerical calculations for 
the orbits throughout their entire evolution Thus our result supports conjecture 2 in (p. 147) 
which states that matter flow lines do not end at a conformal curvature singularity. 

Note that for the type jY models it is only in the SH type V model one has fi = const. ^ 
on the horizon when 1 < 7 < 2 (this follows from Table 8 where we have /i oc e~^-^^^^^+^/^""'''^^). 
Collins has shown that only stiff perfect fluids (7 = 2) allow // = on the horizon |15] in the SH case 
(this implies that there cannot be any fluid flow over the horizon). This lead him to the conclusion 
that this behaviour probably was associated with a rather esoteric equation of state. However, for 
the self-similar models presently considered one always has /i = or /i = cxd at the horizon. When 
the SH type V models are seen as a special case of the present class of models it is rather the case 
IJ, = const. 7^ which is exceptional. Moreover, the existence of whimper singularities is intimately 
connected with // = const. 7^ (see ||T6|| p. 136). Hence we believe that no whimper singularities 
occur in the present self-similar models (a proof would require a more subtle comparison of the 
decrease in and the increase associated with the growth of the affine parameter; see |l^). The 
present analysis (exemplified by the result ^ oc e~^-^(^"*"^+)/(""^-^)) thus reinforce the impression that 
whimper singularities are highly special. 

A Properties of the fluid congruence and Petrov type conditions 

The fluid expansion is given by 

{vv + (1 - v^) {9 - (2a + 3f)Biv)} = 

(35) 

|^ + (l-t;2)(l-(2a + 3/)^t;)|^. 



(1 - 7;2)3/2 



(1 - 7;2)3/2 
The fluid shear scalar is 



^(1 _ ^;2)3/2 



VV 



e 



^(l_^;2)3/2 



^ {l-v^){T.+ -aAv)^ 



{\ — v"^^ (cr+ — aBiv] 



V 



(36) 
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The fluid acceleration scalar is 



3(1 _ 1,2)3/2 



3v + {l- v^) [Ov - 2a+v - 3fBi] 



I {v' + (1 - v'^){v - - 3fA)} 9\ 



3(1 _ ^;2)3/2 

The divergence of the acceleration vector is given by 

p2/x 



+ (1 - v'^){5ve - 4va+ - Sfv'^Bi - 3{2a + 3f)Bi)v 

- (1 - v'^){3fvBi + (2(j+ -e-e^ + 2ea+)v^ - 2Bi{a + f){2va+ - 3fBi))) 



Q.(l-.2)." + i(l + 3.2) {v'Y 



3{l-v'^f 

- ^ (l - -u^) (3Ai(2a + 3/ + fv"^) + - hv) v' 

+ i(l + q){l- v^) ((l - v"^) {3JvAi + 2t;2s+ - v) - vv') 

- (l - v^) {2Ai{a + f)-v) {3fAi + 2v^+ - v) +^ (3fvA'^ + 2v^T,'^)^ 9 
The Weyl scalar is 

- ' -^2 /)_ I o_2 



3^/3 



-3a+ + 3fcS| - 6'cr+ + 2cr1 



26^/" 

3^3 



{S'+ + (2S+ + g)S+ + 3kK] 9' 



(37) 



(38) 



(39) 



Note that the magnetic part of the Weyl tensor is identically zero for all models. The space-time 
is of Petrov type if 

(40) 



- 3(7+ - a+9 + 2(7+ + 3kBl = , 



which can be written as 
otherwise it is of type D. 



(S'+ + (2S+ + g)S+ + 3kK)9'^ = , 



(41) 



B Relation to the fluid approach 

The fluid approach, where one uses a coordinate system adapted to the fluid velocity, is more 
common than the "homothetic" approach used in this article. The line element, of the presently 
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considered LRS models, takes the following form in the fluid approach Q: 

ds^ = -e^WdT^ + e^(^)dX2 + y\\)X^' (dy^ + k'' sin {Vky)dz^) , (42) 



where A = X/T and 6 = {a + f) / f . The relation ak = thus takes the form (6 — l)k = 0. The 
fluid velocity, u", is given by e~*/^(l, 0, 0, 0). 
The following transformation: 

X = e-/(^-^©) , A = X/T = e^* , (43) 

where 

_ = — (44) 

leads to the "homothetic" adapted SSS line element 

ds"^ = e-2^^ds2 = [-N^df + Di^dx^ + Da^e-^'^^ (dy^ + sin {Vky)dz^)'\ , (45) 

where 

(46) 

Dj = /2e-2/^®(e^ - e*-2/*) , = e-^^'^^^Y^ . 
The fluid velocity v is given by f ^ = 
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